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Abstract 

Loop corrections to observables in slow-roll inflation are found to diverge no worse than powers 
of the log of the scale factor, extending Weinberg's theorem to quasi-single field inflation models. 
Demanding perturbation theory be valid during primordial inflation leads to constraints on the 
effective lagrangian. This leads to some interesting constraints and coincidences on the landscape 
of inflationary vacua. 
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I. INTRODUCTION 



In recent years there has been much discussion in the literature about quantum effects of 
long wavelength modes in de Sitter, or slow-roll inflationary backgrounds |1H17|. Depending 
on the authors, these contributions are negligible, infinite, or somewhere in between. A 
clear understanding of these issues is therefore important in light of the experimentally 
verified predictions of the semiclassical inflation theory. Essential to these predictions is 
the assumption that the dominant contributions to density perturbations are infrared finite, 
mode- by-mode. 

In previous work [isl, we emphasized the importance of physical constraints on the choice 
of initial state and explained how this leads to a theoretical uncertainty in the predictions 
for the observations of a local observer. For example, in a global de Sitter spacetime, 
perturbation theory in massive scalar field theory around the Bunch-Davies vacuum appears 
convergent. Nevertheless it is difficult to explicitly introduce an infrared cutoff, and then 
remove it maintaining the symmetries. Depending on one's choice of spacelike slices, such a 
procedure may be necessary. Moreover once massless fields are included (even the graviton) 
the procedure of adopting an infrared cutoff appears to fail, and it seems likely the global 
spacetime is unstable. 

On the other hand, for realistic applications to cosmology we are more interested in a 
local patch of quasi-de Sitter spacetime that expands to our observable universe. In this 
scenario a comoving infrared cutoff is the simplest accurate model, and most of the questions 

nn 

of principle for global de Sitter become irrelevant [18|, 119|. In this context, any sensitivity 
of observables to the infrared cutoff reflects a genuine theoretical uncertainty in predictions, 
originating from the lack of a precisely controlled initial state. Such quantum corrections 



were explored in 



18|. 



In the present work our goal is to extend these results to slow-roll inflation, allowing for 
the nontrivial time dependence of the Hubble parameter. For models built usin g sc alars with 
minimal kinetic termsMt is straightforward to combine the physical setup of 18| with the 
results of Weinberg l|, |2| for this class of models to see that observables at most diverge as a 
power of a logarithm of the scale factor. While this presents serious problems for the global 
stability of de Sitter spacetime, the infrared quantum corrections are tiny for primordial 
slow-roll inflation with realistic parameters. 
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Xue, Gao and Brandenberger 



20| have proposed a related scalar model with non-minimal 



kinetic terms that evade this conclusion. They find large infrared quantum corrections 
produce strict bounds on the scalar couplings arising from convergent perturbation theory. 
In the present work we extend the results of Weinberg to this class of nonminimal kinetic 
term models, and confirm that at most powers of the logarithm of the scale factor appear in 
observables. We then re-examine bounds on the couplings by requiring a good perturbative 
expansion, and find that running of the scalar mass parameters tends to produce an even 
larger effect than that of the infrared modes, with somewhat less strict bounds emerging 



than found in 



20] . We also show that the infrared corrections in slow roll are bounded 



above by the corresponding corrections in pure de Sitter spacetime, as one would intuitively 
expect. 

It is interesting to note that these bounds arising from quantum consistency are not far 
off the kinds of bounds that emerge from tree-level slow-roll considerations, combined with 
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221. We 



matching the scalar potential to the magnitude of observed density fluctuations 
argue this coincidence may be explained using statistics on a landscape of vacuum states. 
Thus the saturation of the perturbative bound on the landscape (at least within this class 
of models) may be regarded a postdiction of the observed density fluctuations. We conclude 
with a brief discussion of how the late-time instability of a de Sitter region is compatible 



with the embedding of a de Sitter region in a unitary model for quantum gravity 
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2J| and 



low the instability timescale that emerges solves the Boltzmann brain paradox of cosmology 



II. IN-IN FORMALISM WITH IR CUTOFF 



We consider slow-roll inflation, with an infrared cutoff imposed as in jl8| . To obtain a 
tractable model of slow-roll inflation we consider a quasi-single fleld inflaton model with two 
scalars: a slowly rolling inflaton (yj) and a spectator fleld (cr), as already considered in 26 1 



and [20]. Using polar coordinates in fleld space, the inflaton and spectator correspond to 
the tangential and radial directions respectively, and the curvature of the inflaton trajectory 
leads to a minimal coupling between the flelds. The action governing the system is 30| 



S 



(1) 



3 



where i? is a constant and the potential V will be constrained in such a way that slow-roll 
conditions for the inflaton are satisfied. We will work in the spatially flat gauge, in which 
the metric is given by 

ds^ = -dt^ + aitfdx^, (2) 

and the scalar metric perturbation has been incorporated in the perturbation of the inflaton 
field if. It will also be convenient to use the collective notation 



^ = • (3) 




Background solution: Perturbing the fields via $ = $0 + and then minimizing the 
action, we find the field equations governing the background solution. Taking the background 
solution to be spatially homogeneous, the background equations are 

<^o + 3i/0o + = 0, (4) 
ao + Sifcro + \/; - ^ = 0, (5) 

where = ^ is the Hubble parameter and = d^V etc. Equation (jl]) places constraints on 
the potential for the field ipo to undergo slow-roll. For the spectator field we pick a constant 
solution, and without loss of generality we can choose ctq = 0. Equation ([5]) then relates the 
steepness of the potential in the radial direction to the speed of rolling infiaton by 

K = f - mr- (6) 

Note the slow roll parameter is non-vanishing for non-zero A 



(7) 



A. The free action 

We wish to use in-in formalism to compute two-point correlators of the form 

where the mode functions 5$ are determined by the free part of the Hamiltonian, and 
the interaction part is taken into account perturbatively as in ([8]). Thus we need the free 



Hamiltonian, which is defined as the part quadratic in perturbations . Expanding the 

action around the background solution $o up to second order yields 



5'free = So + - I d X^f^ 



5'o + - / A; dt 



(9) 



(10) 



where Sq = f d'^Xy/—g[^ipQ — V{ipQ,aQ)] contains the zeroth-order terms, and we have 
switched to momentum space via 



d^k 
(27r)i 



Ak-x 



(11) 



Note that in Fourier space we drop the 5 in front of the perturbation. 
The field equations: From (fTOj) one can derive the field equations 



+ VI - aM = 2A0, 



(12) 
(13) 



We take the infiaton to be massless (V^',^ = 0), and denote by vn? = V'^^ — A^ the effective 
'mass' of the spectator field a. We also take the potential to be of the form V = V{(p) + V{(x) 



to 
of 



eading order, implying V^^ = 0. Both of these constraints are consistent with the analysis 
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26| 



In order to have a solvable system, from now on we will also take the infiaton to roll at 
a constant speed, so A = 0. Then tells us that the potential has to be chosen such that 



the slope is constant along the trajectory. This departs from the analysis of 
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26l. who 



make no such assumption. We choose to set this constraint, because a central tenet of this 

one has to treat A analytically in an 



article is that in order to compute correlators of type (E 
exact manner, as opposed to perturbatively. In 20|, |26| the cross-term Xa(f is bundled into 
the interaction Hamiltonian, whereas we treat it as a part of iffrcc? and restrict to constant 
A in order to be able to explicitly solve the field equations. 

Finally, at the level of the field equations we will work in an 'instantaneously de Sitter' 
approximation, in which the scale factor is given by a{t) = exp[Ht), with a constant H. 



This approximation is valid over time scales 



HhAt < if ^ At < ^' 



47ri?2A2 

which is quite sufficient for our purposes. In particular, it can contain the regime where 
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26| 



effects nonperturbative in A At become important. Such effects are dropped in 
where a perturbative expansion is A is considered. 

At this point it is also convenient to switch to conformal time, defined by 

, , , f dt 1 , . 

dt = a(t)dT, =^ r = —— = (14) 
^ ^ J a{t) Ha ^ ' 

Incorporating the constraints and approximations the field equations in conformal time 
become 

< - + + -WW^) = -^^1' (16) 

where ' = dr- 



B. The free field solution 



We will now solve the field equations fll5fl6p perturbatively in k using the Green's function 
method. Using the expansion 

oo 

s = j2e%, (17) 



i=0 

we can write the field equations in matrix notation as 

= -k'^,_„ with L^\ ^ ^ ' ^ ^ ^' I . (18) 

Note that one should not confuse the mode function $i=o with the values of the background 
fields $0 found earlier. From now on the background values will only appear inside X = (po/R, 
so no confusion should arise. 

We can easily solve $o from L$o = 0, which has power law solutions. One verifies that 
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the general solution is 



$o = XJa(^-J , with a = (0,3,a_,a+), A=\ M, (19) 



3 / . , I f2mV /4A ^ ^ 



where tq is a fixed initial time. We also need the Green's function, defined by 



LG(r,r') = 5(r-r')l, with G=\ ^ . (20) 

\ Ga^ip G(j(T I 

We relegate the computation of G into appendix [Bl here we only present the result 
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Gir,r') = Y,C.Q{r'-r)r' (^)", (21) 

i=l 

where the Ci are constant matrices explicitly given in ( 1B6I) . 

We can now use the Green's function iteratively to solve for higher orders $j. We have 

$,{r) = rdr'G{r,r') ■ $o(r') = {kr^f {-T'^^ " (22) 

Late times: The four independent solutions at late-time take the form 



4X/H \ TO 



(^)" + '*^»''^°" (S)" + I (23) 



where the new coefficients af^ and ^ are functions only of m, A and may be read-off from 
(122|) . The terms • • • denote subleading terms as r — )■ 0. 
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C. Quantization 



We quantize these fields using the mode expansion 

$(x,t) = J rf3ge'«-"<l>(g,t)-a(g) + e-'«-"<l'*(g,t)-a*(g), 

where 

[ai{q), a*{q')] = 6ij6^{q - q') , [^^(g), aj{q')] = , 

and take the vacuum to be the Bunch-Davies vacuum at the start of inflation, annihilated 
by these annihilation operators. In this early time limit, the modes of interest (recall the 
comoving infrared cutoff) are all inside the horizon, and oscillate with time. The quantization 
proceeds in the standard way. 

Later we will need to also estimate the commutator of the fields in the late-time limit, 
where the fields asymptote to the form f|23p . Now the modes of interest are far outside the 
horizon where they decay as real powers of the scale factor f l23|) . At leading order as r — )■ 
(late times), 

where the complex coefficients Cg, Dq, Eg, Fq are fixed by matching to the early time modes 
at horizon crossing. When we compute the commutator [$j(a;, t), $j(x',t')] only cross terms 
between the pairs r°, and r""", r"+ survive, so the commutator falls off as as r — )■ 0. 




D. Late-time limit of observables 

We follow Weinberg jl, 2| when computing the leading late-time terms. As emphasized 
in his work, there are delicate cancellations between terms, which are most easily taken care 
of using the commutator expression 



{Q{t)) = Y.''' "^^^ ^'^^-i • • ■ / [H_ 

N=0 



jih),---[Hi{tN),Qiit)]---]]) 



(24) 



where the subscript / denotes interaction picture operators. Note that rather than taking 
to — !■ — oo as in [1, I2I, we keep it finite, as part of the procedure introduced in [is!] for 
keeping track of the effect of the initial state on observables. Weinberg investigated the 
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leading late-time divergences for massless and massive minimally coupled scalar fields, as 
well as Dirac particles and vector particles. Here we will apply this approach to quasi-single 
field inflation. 

As argued above, the commutator of any pair of elementary fields falls off as or a{t)~^. 
The same will also be true if one considers time or space derivatives of these fields. The 
interaction hamiltonian Hj will be built of products of such fields and derivatives, and 
contain a volume derivative going as a{t)^. As is apparent from (p^ the Hj will always 
appear inside a commutator, so these factors of a{t) will cancel. 

In general observable we may also encounter additional powers of fields that are not inside 
commutators. These introduce at worst constant factors (if only ip{x,t) appears in Q{t)) or 
factors that fall at least as fast as r"-, provided some factors of (y{x,t) appear, or provided 
sufficient inverse powers of a{t) appear in derivatives. 

We conclude then that observables built solely out of products of tf{x,t) do have late 
time divergences, arising from the time integrals in ( l24l) . These lead to divergences as 
powers of logr. These divergences are qualitatively the same as the case of an interacting 
massless minimally coupled scalar considered in [18]. Of course, in order for inflation to 
end, this scalar must acquire a mass, which provides a natural physical cutoff to these time 
integrals. Estimates made in 18| (in the models considered there) show these loop effects are 



negligible for primordial inflation compared to the tree-level contributions, and only become 
significant in genuine asymptotic future de Sitter phases - where they are capable of driving 
an instability. Nevertheless there exists a large class of infrared finite slow-roll observables 
for quasi-single field inflation which contain either derivatives oi ip{x, t) or factors of a{x, t). 



III. LOOP CORRECTIONS AND RENORMALIZATION: INFRARED AND UV 
CONTRIBUTIONS 



The authors of |20| considered a closely related set of questions, and argued much larger 
infrared terms appear due to slow-roll effects. They find divergences that go like powers 
of ( ) ^ /e (with e the slow-roll parameter ([7]), H the late-time Hubble parameter, and 
Hinitiai the Hubble parameter at the start of inflation). Here we will show these conclusions 
change when UV divergences are treated with a physical renormalization prescription. The 
apparent divergences as e — i- disappear, but are replaced by late time divergences involv- 
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Figure 1: Infrared divergent diagram 



ing powers of loga(t). Such divergences match those expected from the de Sitter hmit. 
Demanding that these late-time loop corrections not destroy perturbation theory leads to 
constraints on the interaction potential, which take a similar form to those argued for in 



20|, but differ in the details. 



In section III.l of 20| they consider loop corrections to the scalar curvature perturba- 
tion two-point function (we refer to for details), due to a massless (or sufficiently light) 
entropy perturbation. This can be represented by the a field of the model previously dis- 
cussed, with a higher order ga'^ self-coupling included. The result of their analysis is that 
the leading IR divergences comes from a subdiagram involving the a self-interaction . The 
corrections appears in eqn (56) |^|. Carrying over their result for the subdiagram, and also 
including a physical UV cutoff we obtain 



^^imUaL It ' ^ 2 Of: \ \ n 

Air y \ 1^ initial / J ^iR,phys 'i ■^^ \ \ -'■-'■initial 



+ 9 {^UV,phys ^IR,phys) ; (25) 



where we have estimated the UV and IR divergent terms by breaking the range of inte- 
gration up at a physical intermediate scale, and used the appropriate asymptotic forms of 
the propagator. Note a comoving IR cutoff and a proper UV cutoff is used as in 18|. The 



computation of |20| also uses a comoving IR cutoff, but are less explicit about their choice of 
UV cutoff. A mass counterterm must be chosen to cancel the UV divergence, and impose a 
renormalization condition. This is described in more detail in appendix [XI The result is that 
the would-be Hf^^^^^Je divergence disappears when a physical renormalization prescription 
is imposed for the mass of the a field. 

It is also worth pointing out that even the long wavelength contribution to (125]) is bounded 
from above by the pure de Sitter result, where e = 0. This follows by choosing the scale 
kinitiai = ^iR = C>{Hinitiai) , cLud uotiug for slow-roll with e > 0, the integrand is always 
positive and less than the pure de Sitter answer (e = 0) throughout the range of integration. 
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This point is at odds with the answer obtained in |20|, which may be traced to them dropping 
all but the first term on the right-hand side of (l25l) . 

An important consistency condition is demanding that the perturbative expansion in g 
converge. In analyzing this question, we allow for the mass renormalization of the infiaton 
to incorporate generic effects of new physics near the GUT scale, and still require the pertur- 
bative expansion be valid. The implications for the one-loop renormalization are described 
in appendix [XI Let us examine the physical constraints that emerge from this. If we insert 
;his subdiagram into the full expression for the one loop correction to the two-point function 



18| (simply working in the ga'^ sector of the theory), and ask when perturbation theory is 



valid, we obtain the condition at the end of inflation that 

gH' (iV^ + iV^^) « H\N - log i-Ainro)) ■ 

Now the A^'^ term on the left yields a constraint g ^ ~ 10~^ for massless perturbation 

theory to be valid. The other term requires 

Let us put in some typical values assuming we wish to use the field theory for the infiaton 
from a U V cutoff near the GUT scale (beyond which we expect new physics to set in) , down 
to scales below the Hubble scale. We set HinUiai = lO^'^GeV and Aqut = lO^^GeV. This 
yields 

^7«lO-^ 

I — I 

So we see while these effects are compatible with the bounds of [20] the bounds found 
there are not the whole story, and stronger constraints emerge from the consideration of 
typical UV effects due to renormalization. Another difference with is the powers of A^. 
In both our wor 
elaborated in 



c, and [20| a comoving infrared cutoff is used. The justification for this is 



18] . However the computation of j20[ appears to use estimates for amplitudes 
obtained using the formalism of j8[, who instead use a physical/proper distance infrared 
cutoff. 

Finally we can estimate the two loop contribution coming from the diagram shown in 
figure 2, arising from the coupling of the scalar curvature to the scalar field as described in 
[20]. Following the same type of computation as above, we find a constraint of the form 

gH' (^N' + (%^) + (%^)') « H\N - log(-A.^ro)) . 
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Figure 2: Sub-diagram with two scalar curvature external lines, and the scalar field lines fully 
contracted. 



This then yields the dominant condition 



9N^«1, (26) 



201, though the 



so that g ^ 10^^'^ which numerically is comparable with the bound of 
dominant effect is the short distance renormalization of the scalar mass, rather than large 
distance slow-roll terms. 



IV. COMMENTS ON THE LANDSCAPE 



If we had considered single-field inflation with a potential ga^, a tree-level bound on 
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emerges by matching the observed 6p/p ^ 10^^ with the value predicted by slow-roll 
This yields 5p/ p > Ng^^"^ so that g < 10^^'^. It is interesting to point out that the loop-level 
bound ( l26l) is comparable with this tree-level bound. This coincidence suggests an anthropic 
relation. Namely anthropic/landscape considerations would tend to statistically favor the 
largest value of g compatible with the basic physics of the model. To make g exceed our 
quantum bound requires new physics to appear before the GUT scale, taking us out of this 
class of model. Taking g to saturate the bound, and assuming for the sake of argument that 
one is restricted to working with this family of scalar models, one is then led to a postdiction 
for 6p/p matching observation, jsil 

Finally, it is also worth commenting further on the gravitational version of the late-time 
instability of de Sitter spacetime found in [isl] due to the choice of graviton initial state. It 
was argued there that future eternal de Sitter is actually unstable on a timescale of 10^^^ 
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e-folds. It has already been noted that some kind of late time instability for de Sitter is 



needed for compatibi 



considered in 
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ity with the class of unitary models for quantum de Sitter regions, 



28|. There is was pointed out that this can solve the proliferation of 



Boltzmann brain observers (see 25| for background material). 

In the present context, we find the timescale associated with the production of a Boltz- 
mann brain along the path of some timelike geodesic in an expanding universe to be of order 
its inverse Boltzmann factor e^^^'^ ~ e^^^^ for an observer of order 1 mole of protons, with T 
the temperature of the present cosmological horizon. This timescale is much larger than the 
above instability timescale. The infrared instability of de Sitter spacetime thus has a chance 
to restore our status as typical observers, solving one of the many problems associated with 
doing statistics on a landscape of theory vacua. 
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Appendix A: Renormalization in cosmological spacetimes 

It is helpful to review renormalization in the context of the models considered here, and 



m 
of 



filling in some additional details omitted in the earlier work. The scalar coupling A 
18| will be replaced by g here, to avoid confusion with the discussion in Section 2. The 
loop diagram figure [T] gives rise to the integral of eqn. (9) of [18]. The potentially divergent 
terms arise from the IR and UV ends of the integral, and take the form 




2 



(27r)2ifV4 I 27 2\Hj 2\H) ^ m 



with 7 = m? / (3/f2), and Ty the conformal time of the vertex factor insertion. Here we have 
included the subleading log UV divergent term omitted in [l^, and assumed A^y ^ 

Now let us consider choosing the counter-term so that a physical mass renormalization 
condition is imposed at some scale /i. In keeping with the small mass/early time expansion 
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used in this paper, we expand the IR term in powers of to get 



In this way, we see that with a comoving infrared cutoff, we cannot completely eliminate 
the IR divergence into a mass renormalization, due to the additional logXy dependence, 
which causes problems at very late times. With a proper IR cutoff such a renormalization is 



possible, as discussed in j8|, however as discussed in jl8j such proper IR cutoffs are unphysical 
for spacetimes of cosmological interest (an example being bubble walls moving faster than 
the speed of light). 

Now a mass counter-term produces a shift 

V 

so comparing the UV divergent terms with what we usually have with flat spacetime renor- 
malization we choose to impose the renormalization condition 

Sm' + ^ (a^^ - log ^~H^ log Ain^ = ml,^, . 

Now we wish to impose the renormalization group equation Kuv d^lhysl d^uv = and view 
+ dm? as the bare mass squared ttIq, which implies 

Integrating this equation we find 

ml{Kuv) = ^ (mllogAuv - ^Iv) + c, 

where c is a constant independent of Auy to be fixed by the renormalization condition. 
Substituting we find the physical mass 



m 



phys ^ g^2 

therefore we fix 



c + r-^ {-H"^ log KiR + ml log niphys) 



{-H^ log Air + mlt^y^ logniphys) 



at leading order in g. Thus the loop diagram with mass counterterm gives 
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in this light scalar/early time approximation, and all dependence on the UV and IR cutoffs 
disappears, with the exception of the log(— r^) term which comes from the choice of comoving 
IR cutoff. Note that additional dependence on the IR cutoff appears when the time integrals 
of the in-in formulation are performed, as elaborated in 



Finally, it is useful to reconsider the above computation, assuming instead that new 
physics sets in at some high physical scale Kqut > H. For example, a new field if of mass 
^GUT coupling via ga'^ip'^. In this case we do indeed find a correction of the form 

L(Tv 



showing new physics does indeed lead to a quadratic shift in the mass. This form will be 
useful for estimating the range of the perturbative validity of slow roll theory for energy 
scales approaching the GUT scale. 



Appendix B: Green's Function 

In this appendix we solve equation (l20l) to derive the Green's function of the system. 
Since we know the zeroth order solution (fT9l) . a good ansatz is 

G(r,r') = X^C,e(r' - r)r-^ with = ( ) , (Bl) 

where the coefficients Cij^k satisfy the same relations as 0-^'°^, i.e. 

In order to fix the rest of the coefficients Cij^k we integrate (I20p over the range r G [f'—e, r'+e], 
computing to order C(e°), 

1 = dTLQ{r-r')C,T^'= = - f dr {CuT^'') ~ ^Ir'^- | ^ | "Cfc, 

Jr'~^ Jr'-e ^ - t' \ A oiu - I } 



H 



(B3) 

where summing over repeated indices is implied. The first term arises from c}^0(r' — r) and 
is potentially divergent; we need to demand 

Cfcr'°^ = (B4) 
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for it to vanish. The second term imphes 



2r'"'= I ^ ^ \ ■Ck = l. (B5) 



Expressing C2j,k in terms of Cij^k using f lB2p leaves us with eight unfixed coefficients {cij^k)- 
The remaining constraints (lB4p and fIBSP do not mix Cu^k and Ci2,a:, and hence we are left 
with two groups of four unfixed coefficients, with four constraints for each group. Hence 
solving for Qj^k amounts to inverting 4x4 matrices, and we find the coefficients to be given 
by 



Ck = CkT^-''^ 



m? A 
6{m2+4A2) m2+4A2 \ _ \ 6{m^+4A^ 



^ ' ' 



A 



2A2 2A 



(m2+4A2)V9-4m2-16A2 4m2+16A2-9+3v'9-4m2+16A2 
A(3-V9-4m2-16A2) ^ 

2(m2+4A2)v'9-4m2-16A2 2v/9-4m2-16A2 
2A2 2A 




^ _ I (m2+4A2)^9_4m2_16A2 9-4m2-16A2+3V9-4m2+16A2 i (-na\ 
- I A(3+V 9-4m2-16A2 ) ^ I ' l^^j 

2(m2+4A2)V9-4m2-16A2 2V9-4m2-16A2 

where we set if = 1; it can be restored by scaling m — )■ ^ and X ^. 
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n 

[30] This action differs form the one considered in [26] in that our inflaton has been scaled by 
R for later convenience, ip^s = Rdthem- Note that the kinetic term for d would be —^{R + 
a)'^g^'^ d^Od^O, as expected when 9 is the tangential coordinate and {R + a) is the distance 
from the origin. With our conventions, both scalar fields have units of energy, as does R. 

[31] It is also worth pointing out a similar mechanism may be at work in the Standard Model. With 
the measured mass of the Higgs boson at 125 GeV, the Higgs potential develops an instability 
due to renormalization group flow at an high scale, below the Planck scale [29]. This hints that 
landscape statistics push the parameters of the Higgs potential to a point of marginal stability 
before new physics takes over. 
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